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ABSTRACT 


A remotely controlled wind tunnel model with degrees 
of freedom in roll, pitch, and yaw was designed and con- 
Structed to demonstrate some of the major dynamic stability 
mreecOntrol characteristics of a full scale aircraft. The 
longitudinal characteristics of the model were examined, 
and it was found that the response to a step function input 
mulection of the herrzontal tazl could be» predicted 


accurately. 





TABLE OF CONTENTS 


ee ene DOGG) te = = = = = = = 
Pe NODE CDE SeCRIRITON= = 9- -—=- - - 
Oe NEw ciEaniKEs = =e 

B. PHYSICAL MODEL DESCRIPTION - 

HE USC diem inge a. = f= > = 


2. Wing ---------- 


See HOT LvoOntal and Vertical Tails- 


DN our rORT (BEARING == "= - 


C 

ree CE UNAS MBLY = 9-2 >= °- 
IE ee LUNN Ee MOUNT === = = = 
|e} 


SUMMARY OF MODEL DIMENSIONS- 


Pew we ieCAl MODEL - --- = = = = - 


A. LONGITUDINAL EQUATIONS OF MOTION 


B. LATERAL EQUATIONS OF MOTION- 


ieee tABDEITyY DERIVATIVE RELATIONSHIPS - 


ae NGrTGDINAL STABILITY DERIVATIVES 


ProtennL stab lia iy DERIVATIVES- 


V¥V. EXPERIMENTAL LONGITUDINAL RESPONSE 


A. EQUIPMENT CALIBRATION- - - - 
POESEDEO@CRDIIRE = = = = = - = 


Coco ens tON se e= = - = = = = 


—_ 


Vie RECOMMENDATIONS FOR FURTHER STUDY - 


-_ 


APPENDIX A: LATERAL STABILITY DERIVATIVE EQUATIONS- 


COMPUTER OUTPUT- - - - - - - - - = - - 


a 








COMPUTER PROGRAM - - 
IpomeOr wR ErERENCES = 


INITIAL DISTRIBUTION 


Bod = 





Lis Bee 


TJs ie 


let 2. 


ii 3. 
i= 4 
= 5. 


i - 6. 


i 7. 
iS. 


Vaan 


Ria « 


bis! OF FIGURES 


Medei@tcacmemceyotability Axis System- - - - = -=473 
Model Configuration- - - - - - - - - - - - - - -20 
Model Internal Components Showing the Receiver 

and Horizontal and Vertical Tail Servos- - - - -21 
Motion Sensor Circuitry- - - - - - - - - - + - -22 
Wanineeocetul Oil amem- - = = 2 = = = = = =< = << = = = 25 
Main Support Bearing - - - - - - - - - --+-- -24 
Tail Unit Assembly with Tail Surfaces 

LinS Ga SG= Se A 
Wind Tunnel Mount- - - - - - - - - - = - = - = -26 
Model Mounted in the Wind Tunnel - - - - - - - -2/7 
Poucwmmloneter Galtpration Curves = = = = = = = =47Z 


Control Surface Angle vs. Transmitter Stick 


rete eee ee ee = ee ee a A? 
Angle of Attack vs. Time, wate SUR ett, See. a) 45 
Anele of Attack vs. Time, Vee $8.7 ft/sec - - -44 
Angle of Attack vs. Time, ue = 65.4 ft/sec - - -45 
Angle of Attack vs. Time, Ws. 71.7 ft/sec - - - 46 
Undamped Natural Frequency vs. Tunnel Speed- - -4/7 
Damping Ratio vs. Tunnel Speed - - - - - - - - -47 
Undamped Natural Frequency vs. Friction 

Coefficient- - - -----+-+-+-+--+-+--- + -48 
Damping Ratio vs. Friction Coefficient - - - - -48 








biol OP yor MBOLS 


fitt curve slope 

HomaZoneda tall Jiftt curve slope 
Wing span 

Mean aerodynamic chord of the wing 


Vaniatron Of seollling moment coefficient 
haee pentulbpations 


Lomlairon) Cn holling moment coefficient 
yaw rate perturbations 


Variation of rolling moment coefficient 
Sideslip angle 


Vaoutaevon Of trolling moment coefficient 
Veomtlred! tari perturbation angle 


Variation of rolling moment coefficient 
aileron perturbation angle 


with roll 


with 


with 


with 


With 


Variation of pitching moment coefficient with 


pitch rate perturbations 


Variation of pitching moment coefficient with 


age Ot attack 


Variation of pitching moment coefficient with 


time rate of change in angle of attack 


Variation of pitching moment coefficient with 


eaimeeanele or attack perturbations 


Variation of yawing moment 
roll rate perturbations 


Variation of yawing moment 
yaw rate perturbations 


Variation of yawing moment 
Sideslip angle 


GOcrtricilent With 


coefficient with 


CoOcrtrmetent with 


Variation of yawing moment coerficient with 


Vettveal tail perturbation angle 





> - oe > — 


op SMe oo ae © f) 


= = &> 


Variation of yawing moment coefficient with 
aileron perturbation angle 


Demeresopyenc derivative of a quantity with 
respece tO NnOndimensional time 


C.G. position, fraction of the mean aercdynamic 
chord 


Netamalenornteor the modell, fraction of the mean 
aerodynamic chord 


Nondimensional moment of inertia of the model 
apoE une XxX axis 


Nondimensional moment of inertia of the model 
apout the Y axis 


Nondimensional moment of inertia of the model 
abouteche 2 axis 


Nondimensional product of inertia of the model 
about the Y axis 


Dismoimee trom the C.G. to the horizontal tail 
aerodynamic center 


Dicmanec trom the C.G. to the vertical tail 
aerodynamic center 


Roll rate perturbation 

Pitch rate perturbation 

Vanerate perturbation 

Wing area 

HOtmmZzGmtal: tall area 

Vertical tail area 

ESmuliaiattons of airspeed in the X direction 
Hommuzomtar stall volume ratio, OS aCe s} 


Vem@mrecal tall volume ratio, Cs) 7 (© Ss} 





Greek Symbols 


oO Angle of attack 

a, Perturbation of the horizontal tail angle 
B Sideslip angle 

op Perturbation of the vertical tail angle 

G Damping ratio 

n Petelroatronsarethne elevator angie 

u Ppletuon COet ri cient 

E Pertui mation of the aileron angle 

O Undamped natural frequency 

superscripts 

° Variation Of the quantity with respect to time 


10 


ACKNOWLEDGEMENT 


The author would like to express his thanks to 
Mr. Ronald Ramaker and Mr. Glenn Middleton of the 
Aerospace Department's wood and metal working shops for 
myeir time and effort in the actual construction of this 


model and its related components. 


id 








tea ANT RODUGITON 


Mimtncronlmuoneot the author, in order for a student 
to realize full benefit from any classroom learning process, 
that process should contain, whenever possible, some type 
of functional visual aid or laboratory to supplement and 
Support the in-class lecture periods. The purpose of the 
design and construction of the stability and control wind 
tunnel model described in this thesis was to provide that 
additional support element to the stability and control 
courses offered at the Naval Postgraduate School, Monterey, 
Salrcornia. 

Although it was not possible to incorporate all of the 
degrees of freedom of an actual aircraft into the model, 
the model was designed to demonstrate many of the major 
femamnre characteristics of a full scale aircraft. The 
response of the model differs somewhat from that of an 
actual aircraft due to the fact that it must be supported 
within the wind tunnel test section and prohibited from 
translating in any of the three coordinate functions 
(Fig. I-1). This difference, however, is slight and does 


not detract from the primary purpose of the model. 
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Figure I-l. 


Model Reference Stability Axis System 
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II. MODE DESCRIPTION 


A. GENERAL FEATURES 

Pacwmodelm(iag. 11-1) was désigned to illustrate the 
characteristics of a typical single engine mid-wing aircraft. 
With minor modifications, however, the effects of various 
wing geometries can be studied because the model was designed 
With a removable wing section. Similarly, model tail sizes 
and shapes are variable. By the use of various wing and 
Poe geometries the Stability characteristics of the model 
may be changed. 

Model attitude was controlled by the deflection of 
ailerons, and by horizontal and vertical tail deflections. 
The horizontal and vertical tails were designed as "all 
flying" surfaces so as to eliminate the need for Separate 
elevator and rudder. sections. The control surfaces are 
controlled from outside the tunnel by means of a remote 
control transmitter. A receiver mounted within the model 
receives the signals from the transmitter and, by means of 
mechanical linkage, a servo deflects the appropriate surface 
mer Control (Fig. II-2). 

The resulting model motion was detected by three 
potentiometers which were situated so as to determine the 
emeilar rotation Of the model about the X, Y, and Z 
Stability axes. The potentiometer output was connected 


to a chart recorder which supplied a graphical history 


14 





Pmevumesmodeleresponse. Ihe motion Sensor circuitry is 


Shownoein Fieure L1-3. 


B. PHYSICAL MODEL DESCRIPTION 
1. Fuselage 

The fuselage has an overall length of 29 inches 
and was of circular cross section. It has a maximum out- 
Side diameter of 4-1/2 -inches and a fineness rate of 6.44. 
The fuselage 1s constructed from a block of mahogany which 
was hollowed to give a wall thickness of 1/4 inch. 

Ze Wome, 

The wing 1s constructed of mahogany and has an 
NACA 2415 airfoil section, 6 degrees of dihedral, and no 
geometric twist. The wing has a span of 30 inches, a 
leading edge sweep of 23.5 degrees, and a taper ratio of .3. 
The left and right wing panels are attached to the fuselage 
Ppeemeans Of a carry through section (Fig. II-4) which also 
houses the main support bearing (Fig. II-5). The rectangu- 
Pnoles in the carry through section were cut to reduce 
ne weight of the unit. 

Ailerons, which are 8 inches long and 1 inch wide, 
were cut from the trailing edge of each of the wing panels. 
ime total aileron area is 13.33 per cent of the wing area. 

oe WOneZOncaleandevertreal Taiwls 

Both horizontal and vertical tail surfaces are 
constructed of mahogany and have an NACA 0012 airfoil 


section. The tails are designed as slab surfaces so as to 
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eliminate the necessity for separate elevator and rudder 
sections. The tail surfaces are held stationary on rotating 
axles by means of set screws. The axles are in turn sup- 
ported by the ball bearing tail unit assembly shown in 
Paeure LI-6. 

The horizontal tail has an overall span of 12.8 
inches with a ieading edge sweep of 25.2 degrees anda 
Paper ratio of .24. The horizontal tail area is 26 per 
cent of the wing area. 

The vertical tail height (from the fuselage center- 
line) is 6-1/2 inches, the leading edge sweep is 50 degrees, 
Mmamenice taper ratio is .31. A vertical tail area to wing 


‘mEcaeratio vor .22 15 used. 


C. MAIN SUPPORT BEARING 

The model is supported by the gimbaled arrangement 
shown in Figure II-5. This unit is installed so that the 
Pivot point of the bearing is located exactly 14-1/2 inches 
aft of the nose. This point serves as the origin for the 
model reference axes aS well as the model center of gravity. 
Potentiometers incorporated within this unit detect model 
development in roil and pitch up to angles of 25 degrees. 
These potentiometers are used as voltage dividers, and their 
SmtEput voltage 1s connected to a chart recorder which pro- 


vides a graphical measurement of the model response. 
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Dee TALL UNIT ASSEMBLY 

The assembly shown in Figure II-6 serves to support the 
horizontal and vertical tail shafts which rotate the sur- 
faces to the desired deflection angle. The unit is con- 
structed of aluminum and contains four ball bearings (2 
mem izOntal and 2 vertical) on which the tail axles are 
Supported. The Gntire assembly is attached to the fuselage 
at a point such that the horizontal axle is exactly 12 
inches aft of the model center of gravity. This resulted 
in a horizontal tail volume ratio, Vip Ole 2 aCe 


vertical tail volume ratio, Vy» Oto. 


E. WIND TUNNEL MOUNT 

The model was supported in the wind tunnel test Section 

by means of the mounting assembly shown in Figure II-7. 

This support was constructed of aluminum and was attached 

to the roof of the wind tunnel as shown in Figure II-8. 

It has an overall length of 12 inches of which approximately 
Six inches protruded down into the tunnel flow. This portion 
of the support was shrouded by a symmetrical aerodynamic 
falring . 

The model was connected to the mounting unit through a 
7/16 inch diameter stainless steel tube. The Shaft of the 
main support bearing was inserted into the bottom of the 
tube and tightened in place with two set screws. Collars, 
in conjunction with thrust bearings located within the main 


Miiiewaosorped any lift, drag, and side forces of the 
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model while still allowing it to rotate freely about all 


the coordinate axes. 


The shaft of a potentiometer located on top of the 


mounting unit was inserted into the top of the stainless 


steel tube and was held tight with a set screw. 


This 


potentiometer sensed yaw angle of the model, and a yaw 


meopper attached to the steel tube limited yaw to 30 degrees. 


F. SUMMARY OF MODEL DIMENSIONS 
Wing: 


Area (including portion covered 
by fuselage) 


opan 

Cond -root | 
aaa 
-MAC 


Mepect Ratio 
Momnotl vsec t10n 
Lere ane Edge Sweep 
Dihedral Angle 
Homezontal Tail: 


Area (including portion covered 
by fuselage) 


opan 

Cnonma-1Too0t 
-tip 

Aspect Ratio 


Peto le occ t1on 
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Leading Edge Sweep 

Incidence (variable) 

Tail Volume Ratio 
Mertacal Tait: 


Area (including portion covered 
by fuselage) 


Height (from fuselage centerline) 
Chend=noo t 
=) IE 
Airfoil Section 
Leading Edge Sweep 
heme OluUMNe) Ranta © 
en Die SLeM ere + 
Total Length 
Maximum Cross Sectional Area 


Fineness Ratio 
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III. ANALYTICAL MODEL 


The model was designed, and its motion predicted by 
Heans Ot linear CGuations of motion. From these equations 
it was possible to predict model response to various control 
surface inputs. Also, once these equations had been estab- 
lished, it was possible to theoretically vary model charac- 
teristics and obtain information on how this new configuration 


would respond. 


A. LONGITUDINAL EQUATIONS OF MOTION 

The nondimensional longitudinal equations of motion 
were obtained from Ref. 1. These equations were developed 
using the small disturbance theory for an aircraft having 
four degrees of longitudinal freedom. The model, however, 
was suspended in the wind tunnel so that translaticns along 
the X and Z axes were denied, thereby allowing the X and Z 
force equations to be neglected entirely. The elevator 
hinge moment equation was also neglected, as control dis- 
placement rather than control force was considered as an 
input. 

Due to the assumptions of smali angles and constant 
airstream velocity in the tunnel, perturbations of the 
meoleocity in the X direction, U terms, were zero. AI1so, 
the model was situated within the tunnel so that the flight 


path angle, 8, and the angle of attack, a, were equal. 
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Die TOME ncCuctE tide the horizontal tail was a slab 
surface and tnere was no Separate elevator section, the 
quantity n which was originally defined as elevator per- 
turbation angle was replaced by tail angle of attack Cy 
It was further assumed that a terms were negligible. 

Applying the above assumptions to the initial equations 


of motion yielded the following longitudinal equations of 


motion for the model. 


ipa - (one to C oS ce a - C a, = 0 III-1 


q= 6 ey. 


These equations, however, do not take into account the 
Friction inherent in the main support bearing on which the 
model pivots. Therefore, they must be further modified to 
account for this frictional retarding moment. 

The frictional moment present at tne pivot is a function 
weemodel Supporting force. More accurately: 


liber eee, Sue) Inte 3 


where: Wy = @eerricrent Of friction {1n.) 


model weight (l1bs.) 


=, 


= 
I 


mogel lift (1bs.) 


EUttimewtnis equation in nondimensional format results 


mn the following expression: 





C _ 2 “ac Mu Le] C OL 
io (i 
fe rope Wi mosve C O 


IIt-4 


where: op = air density (slugs/£t”) 


V 


velocity (ft/sec) 


Maeetirst term on the right hand side of equation III-4 is 
the constant break-out moment of the bearing for any given 
airspeed; therefore, the perturbation of the frictional 


moment coefficient during motion was given by: 
i 1 7 
AC = O a as 
m 


and the revised model longitudinal equations of motion are: 


Qe es. 


u 
es ee 2 1 % 
ipa - (Cn 7 C )a- (Cc - macy cer Cs a 


q O Cc OL mn 


The basic rules for the formulation of state equations 
from ordinary differential equations [Ref. 2] were followed 
inmeerder to obtain the state variable representation of 


the model longitudinal equations of motion: 
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B. LATERAL EQUATIONS OF MOTION 

As with the longitudinal case, the lateral nondimensional 
equations of motion were obtained from Ref. 1 and then modi- 
fied to fit the explicit conditions of the model suspended 
within the wind tunnel. To begin with, the Y force equation 
was disregarded due to the fact that the model was restricted 
fiom translating along that axis. Also, the rudder and 
aileron hinge moment equations were neglected since control 
displacement rather than control force constituted the input. 
The term ~ in the original set of equations denoted rudder 
omele perturbation; however, this was replaced by Sy since 


miresvertical fin was designed and constructed as a slab 


Surface. It was further assumed that all E and é.. terms 


F 
here zero. 
The model was suspended within the tunnel so that the 


rate of change of the bank angle, b, and roll pate. D, were 


si) 





equal. Also, the rate of change of the sideslip angle, B 
was the negative of the yaw rate, r. 
Applying the above conditions to the initial equations 


of motion resulted in the following model lateral equations 


Oe mMOtLOn. 


ate > + (1, - C, )p - (i,D + ad - C, EF - C 6. = 0 III-9 


1, E : 1g. F 
fmee = (i-D + C.. )p + GUD Scie \r - CS SCR vege oti) 
B p r E Sy 
$6 =p  III-11 
g =-r Wee 


However, as before, these equations do not reflect the 
frictional moments present and must be modified to do so. 
Model side force is the major parameter contributing to the 
retarding frictional moments during model lateral motion. 


Therefore: 


Lerict. ~ 41" go oe 


and 


Neto lene Tii-14 


~JiGretiomal conponent: Of rolling moment 
(in-lbs.) 


where: ee 


= frictional component of yawing moment 


No. 
PMCS eS 
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Hy > = COCmeNeTeMes Of L£riceson (in. } 
b | 


ve = side force (lbs.) 


The frictional coefficients Wy and Lo Should be equal, but 
the coefficient for the yaw axis has a different value due 
to the contribution of the yaw potentiometer. 

Putting these equations in nondimensional format results 
mamtne following expressions for the respective frictional 


moment perturbations. 


C eltcne I 
A : C TI-15 
and 
Vo 
AC = 5 C B ITII-16 


where: b = wing span 


Incorporating these contributions to the rolling and 
yawing moments into the previous set of equations results 
in the following corrected set of model lateral equations of 


motion. 


u 
1 ; A | A 
- (Cy + eye GDN CP 2) (CLD ONG Ds 


Gee y. Guo >. 0 Tey 
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Ng 6 r 
- * - os b= 0 Telele=ain 

F 
$ = p I Tieko 
See ee III-20 


Putting these equations in standard state variable 


format produced the following matrix system. 
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Lee ome DY Oh VAT VE. RELATIONSHIPS 


In order for the equations of motion to be of any use 
in the design of the model, the individual stability deriva- 
tives were programmed and expressed in terms of the model 
physical characteristics. The following sections present 
the relationships used to compute the required stability 


derivatives. 


PeeeeLONGITUDINAL STABILITY DERIVATIVES 
The values of the stability derivatives involved in the 
longitudinal system of equations were determined using the 


relationships presented in Ref. 1. These are: 


C = -a(h_ - h) IV-1 
m n 

OL 

1, 

C - ~2a,(—)Vy, IV-2 

q e 
an ae G& IV-3 
m° il (aves 

q 

“mn, = a, Vu IvV-4 
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From the above equations it is evident that the longi- 
tudinal characteristics of the model are very dependent 
upon the tail length (1,) and the tail size (Si), dda. t nat 
Dy Waryine the tail area, the longitudinal stability charac- 
teristics of the model may vary significantly. For the 


model under study in this paper: 


ee 
I 


=A slatat Come l4. 979 
q 


Ores oul C = -2.9448 


om 
UI 


bee LATERAL STABILITY DERIVATIVES 

Unlike the longitudinal case, the expressions for the 
lateral stability derivatives are very complex and unwieldy 
to use and program. The equations of Ref. 3 were used to 
estimate these derivatives and are given in Appendix A. 
For this particular model, the computer program solution 


had the following values: 


Cy = -,1254 C. = .07379 
B B 

C = -,.4381 C = -,00251 
a n 
p p 

Cy = .0437 C = -,15384 
ie ny. 


oy 
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V. EXPERIMENTAL LONGITUDINAL RESPONSE 


A. EQUIPMENT CALIBRATION 

With a 4-volt potential across the sensor system, the 
Pech and roll potentiometers were calibrated using a 
clinometer, and the yaw potentiometer was calibrated using 
a protractor with the yaw stopper as a pointer. The result- 
ing calibration curves are presented in this section as 
hoeoure VYV-1. The relationships between control surface 
deflections and transmitter stick angle were also measured 


and are presented in Figure V-2. 


Pree lest PROCEDURE 

The sensor system was energized, after which the wind 
tunnel was activated and set at the desired airspeed. The 
model was then trimmed out at a slight nose-up attitude. 
The horizontal tail surface was deflected and the model 
response to a step function input was observed and recorded. 

Losterinicwat airspeeds Of 50.7, 56.74.6524, and 71.7 
ft/sec and tail angle changes of 3, 4, 5, and 6 degrees 
(positive tail angle is with the leading edge up) were 
accomplished. The graphical results of the model response 
to a tail deflection of 6 degrees are presented as typical 
mie rigures V-3 through V-6. 

From these results, it was possible to determine the 
values of the damping ratio, g, and the undamped natural 


frequency, Wa? for the model at the different airspeeds. 


Se) 





beapires V-7 and V-&8 clearly illustrate that as the tunnel 
wpeedmls imc¢measeda, the value of om Temas Cairnky constant, 
Wide the value of ¢ decreases. This results in a gradual 
decline of the damping factor 25W,, as the airspeed is 
increased. This decline in the damping term is a conse- 
quence of the increased wing lift at the higher airspeeds 
which reduces the pivot force resulting in less bearing 
felctlon. 

Equation III-6 was used to obtain the characteristic 
equation for the model and the following expression for 
w and z were obtained as functions of the friction 


n 
coefficient, Wy: 





Figures V-9 and V-10 illustrate these relationships for 
this model assuming constant Wy values. 

The experimental] results were used to determine the 
maknown bearine friction coefficient, Hy - It was found 
Mietea value of 1.75 was required to match the predicted 


and experimental values of Wo. Figure V-8 snows the 
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resulting comparison values for the damping ratio. It is 
evident that as the tunnel Speed increases the experimental 
and theoretical values show fair agreement. 

From the various test runs completed it was found that 
tunnel speeds around 65 ft/sec exhibited the best charac- 
teristics for both visual observation and graphical data 
acquisition of model response. At velocities greater than 
75 ft/sec considerable model buffeting occurs which tends 
to overshadow any data obtained. Conversely, at tunnel 
speeds under 50 ft/sec the model response is greatly in- 
habtted by bearing friction, and the resulting heavily 


damped motion rarely has more than one oscillation. 


eee CONCLUSIONS 

The model studied has demonstrated that the basic 
concept of a remotely controlled stability and control wind 
tunnel test model is a valid one. It has been shown that 
bometne longitudinal case, the damping and frequency of the 
model may be predicted quite accurately using linear 


equations of motion. 


al 
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Figure V-1. Potentiometer Calibration Curves 
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Figure V-3. Experimental Angle of Attack vs. Time 
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Figure V-4. Experimental Angle of Attack vs. Time 
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Passure V-5. Experimental Angle of Attack vs. Time 
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Figure V-6. Experimental Angle of Attack vs. Time 


via = 71.7 ft/sec, Ao, = 6 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 


Mie tO Insuitici1ent time, an extensive study of the 
lateral characteristics of the model was not possible. 
However, it must be pointed out that the model exhibited 
very unstable handling qualities in roll. In fact, this 
instability was detrimental to the point where the model 
had to be restrained from rolling in order to obtain the 
degree of accuracy required to satisfactorily determine the 
model longitudinal characteristics. Further study in this 
area of the model's response may reveal the reason for, as 
well as a solution for, this unpredicted instability. The 
model lateral response should be carefully examined and 
compared with the theoretical values obtained from the 
computer program. 

The fact that there presently exist difficulties in 
the roll characteristics of the model should not immediately 
bring a conclusion that the basic concept of a three- 
degree-of-freedom stability and control model is impractical. 
Rather, with continued testing and refinement it is believed 
that a tool such as this will become an invaluable asset 
to students actively engaged in the study of aircraft 
Stability and control. 

Additionally, it is suggested that the friction 
coefficient of the main support bearing and yaw potentiometer 


be experimentally calibrated so that the model response may 
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be correlated more accurately. In fact, a total redesign 
of the main support bearing for less friction would greatly 


simplify and increase the accuracy of the predicted response. 
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APPENDIX A 


CAE Reies baie ty DERIVATIVE EQUATIONS 


The following equations for the lateral stability 
derivatives were obtained from Ref. 3, and, unless other- 
wise stated, any specific figures mentioned in this appendix 


mre from that reference. 


Pee vARTATION OF SIDE FORCE COEFFICIENT WITH SIDESLIP ANGLE - 


é 
YB 


The wing contribution to this derivative may be estimated 


from: 


C = som. |i) es as A-2 


Meme | 1S the wing goemetric dihedral in degrees. 


The body contribution may be found by using: 


a wing-body interference factor obtained 
oy Jaa yal 


It 


Where:  K. 
1 


S. = the cross sectional area of the fuselage at 
the point along the body where the flow 
ceases to be potential. 


HiCElbeGataereLerence area. 


Cp) 
Il 
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The vertical tail contribution may be estimated uSing: 


5 
— do BAN - 
°y k cL (1 + dB) Ny xX A-4 
V V 
where: k = an empirical factor defined in Fig. 7.3 
Cy = vertical tail 1ift curve slope 


Se Z 
do 7 V W 
Clie aa? Nyy ae eos F500 iP aaa + 4 APs Pee OOS ex 


Y 
il 


isn ermalic 2) ll wyers tou) lee parse! 


V 
> = reference area 
A = aspect ratio 
Z = vertical distance from the wing root quarter 
chord to the fuselage centerline (positive 
for wing above centerline) 
d = fuselage diameter at wing intersection 


Pee vARTATION OF ROLLING MOMENT COEFFICIENT WITH SIDESLIP 
Pavia = 9G 


The wing-body contribution may be determined from: 
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“al 


where: 


BWR 


“1, “1, 





E 
ie 
Some int abs K. + (a —) } + T[5— Kk 
Lip Cy Ney 2 M, f£ A It Mp. 
AC 
& By co GG. peas, A-6 
| ays 
B Ww 
Cy = wing-body lift coefficient 
WB 
eal 
—), /2 = wing sweep contribution obtained from 
Dae Ralee we 7 eelel 
Ky = the compressibility correction to sweep 
A ete Oil yee Day ee 
Ke = the fuselage correction factor from Fig. 7.13 


1 
(rc), = contribution from aspect ratio obtained 
L fae Omen dete 17 2.14 
r = dihedral angle 
oF 
—* = dihedral effect of the wing from Fig. 7.15 
Kay = compressibility correction to dihedral 


r romero, 7216 


mr d, 2 

27 = U0) (5) 

b = wing span 

d = diameter of Sonenne at wing intersection 


mL Zk 
a oe es SneS, a3 ( 


ayes 





ZT was previously defined for equation A-4. 


mie horizontal tail contribution is given as: 


ay i PH OH A-7 
H HB 
where: Su = horizontal tail area 
bn Homizontal staal span 
Cy 1s obtained by solving equation A-6 using the 


Bue HOt~MZontal tar) enareaeteristies vice the wing. 
Mme vertical tail contribution is given by: 


Z 1 : 
: (Gee Osmo — 6 1VaSin. dd) 
G C a 
Yg b 


where: oo was previously determined using equation A-4 
By 
a = angle of attack of the wing 


Z.. = absolute distance from fuselage centerline 
COUVeLtlcalstalilvacrodynamic center (a.c.) 


iw IpsOllte distance from model c.g. to vertical 
ail. 


Sree oaRITATION OF YAWING MOMENT COEFFICIENT WITH SIDESLIP 


ANGLE = C 
Ne 


The body contribution to equation A-9 is given by: 
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Cp) 


one 
C. a ces eS Ry Kp EE ae A-10 
By ] 
roe Te ; Ky = an empirical factor for body and wing-body 
CpeeGus trom bre; 7.19. 
Kp = Reynolds Number factor for the fuselage 
It i Ole erm? ec.0 
Sp = body side area 
S 
1, = overall fuselage length 
Mie vertical tail contribution is: 
C ee (ty cos a + “vy sin ao) A-11 
"8 ”8 
V V 
where all the elements have been previously defined. 
eee VARIATION OF ROLLING MOMENT COEFFICIENT WITH ROLL 
RATE - Cc, 
p 
Cy = C, e Cc, 1p Cy joa i 
P PwB Pi Py 
ives wing-body contribution to the above equation 1s: 
BC, 
e k 
e ~ C = (|) B A-13 


1 i 
Pw Pw 


BC 
Where: ( ) = roll damping parameter from Fig. 8.1 
p p 


ae 








k = the ratio of the average wing section lift 


curve slope Cy to 27 


Mie horizontal tail contribution is: 


ea ~ 8 oT dy So (5) 
Py Pp 


where: (C, dy 1s obtained by solving equation A-13 for the 
p Hora zZOnmtad.” taal 


The vertical contribution is estimated using: 


a Ly, Z 
C = 2 Ge C A-15 


iL 


where all the variables have been previously defined in 
this appendix. 
E. VARIATION OF YAWING MOMENT COEFFICIENT WITH ROLL 


male - C 
n 
p 


The wing contribution is estimated using the following 


equation. 








n 
— ~- = ~ =- ~ 
C Cy tan a [ C, tan «a Cm) r CJ A-17 
Pw Pw P eae 
M 
where: 
Cy PeetOund £50M equation A- 15 
Pw 
® = wing angle of attack 
Cy = wing lift coefficient 
C il Z C 
n : a AB+>(ABtcosA. ,,) tan Maya 2) i 
C = +4cos All ee - 
ie Cy 0 ey aa Doyles *cosA.,4)tan Naya ic Cy 0 
M M=0 
feee= aspect ratio 
re } 
pee \/ |-M cos ey 
M = mach number 
2 
_ tandA tan A 
C A+6(AtcosA_,,) (= ome 2 ee 
n 1 C/idare A Ane 
(c*) ee ic osiiis A-19 
I, C,=0 e/a 
M=0 
The vertical tail contribution may be found from: 
7 . Z.,cos Gas 1 sin 7 oe = 
Cc aes (1 cos Gat Z S10 ee) ae ) “y, a 
Oey) V 


where all the elements of A-20 have been previously defined. 


Sa 





F., VARIATION OF ROLLING MOMENT COEFFICIENT WITH YAW 


KATE = Cy 
Yr 


The wing contribution is given by: 











C AG 
= eS r 2 
oT CT (ee a0 * OFT) a 72 
Ty L L 
M 
where: Cy = Wing lift coefficient 
2 AB+2cosA_,, tan“A 
1 2B(AB+2cosA ABt+4cosfA 8 1 A-23 
(ay _ ine a SO 
i “L=0 At2cosi_,, tan“A CL “¢, =0 
M 1 ei 4 c/4 
At+4cosA 8 M=0 
Gy 
3 
—) is obtained from Fig. 9.1 
L C,=0 
M=0 
“aa Asink. 4 
fF 12 At4cosh. yy 


T = wing dihedral angle in radians. 


ine vertical: tail contribution is found from: 
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Tk : : 
Fi. re: (1,cos a + Z sin a)(Z cos a - 1 sin a) Y A-24 
J V 
where all the variables have been previously defined. 


G. VARIATION OF YAWING MOMENT COEFFICIENT WITH YAW 


RATE - C 
ny 


G 6 
C = ans, co 2) C A-26 
n,. C 2 i Cy Da 
W L O 
where: 
C 
My 
a) eS SEONG ee mom Lo Oo 4c 
oF 
C 
ae 
cree 1S) OMMC moO bao. DSi. 
Do 
Cy is the model zero-lift drag coefficient. 
O 


Tne vertical tail contribution may be estimated from: 


= 2 - 1 2 Pans 
Co 5 (1, cos a + 2 sin a.) Cy A-27 


where all the variables have been previously determined. 
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H. VARIATION OF ROLLING MOMENT COEFFICIENT WITH AILERON 


Pier! ON = C 


1, 


The procedure for determining this coefficient is a 
rather complex step-by-step process highly dependent upon 
fmaemecxact type Of aileron employed. Therefore, it is sug- 
mested that Ref. 3 be consulted for the determination of 


mrs derivative. 


I. VARIATION OF YAWING MOMENT COEFFICIENT WITH AILERON 
VE becrion = - os 


E 
C = Gao A- 28 
Ne L 1, 
where: K —Boilecnprer teak actor Obtained £rom Fie. 11.3 
Cy “=oleeMmodelus teddy state ditt coet ficient 
C, = the derivative determined in section H 


J. VARIATION OF ROLLING MOMENT COEFFICIENT WITH RUDDER 


DEPLECTION = Cy 
OF 


Reference 1] was used to determine the value of this 


derivative. The following equation was used. 





where: Cy Veruredime tall lift curve siope 


Wekeledietai | area 


Up 
il 
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K. VARIATION OF YAWING MOMENT COEFFICIENT WITH RUDDER 


DER EEGEION = Co 


oz 


Reference 1 was used to determine the value of this 


derivative also. The following equation was used: 


minere all the variables are known. 
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Fok kok i kok lot dak gota iotk fea kde tek dak kkk ee ak aa ak tk keke ok akc ake ok fe akc ake ak ofc ake ak ofc ake ake of of 
Fe ek Xe fe fe 2k ke ak of ok 2 eee PORE ea eeeueanenaes cua e 


xK XK KK 
XK XE aK Xs 
KX MODEL CHARACTERISTICS %K AK 
mK XK x 
% XK KE 3k 
ME He OK Me ake Oe KK kK NS OK Ne Ne i NK eC AK OE OK AK EO NE OI EK AK OK OK OK Ss aad SD) MK XK OK FS 
EK EK KS KK ON OK KK OK OK OK KK eX KE OK KE KX NEE oie OE KE KK Oe Eg ie Xk cK OS aK ok ke oie se ie ook aX 
TOA MODE SPAR AMETERS = 
tciMeelet COEFFICTENT (CL) = 0.0 
MODEL ANGLE OF ATTACK (DEG.) = -5.0982669)0-O0l1 


TAIL ANGLE OF ATTACK (DEG) 2647659510 00 


TOTAL LENGTH CIN) = 2.9000000D Ol 
TOTAL WEIGHT (LBS) = 4.6997845D GO 
Sl Cire ey STATIC MARGIN: = 222599784D-01 
LIFT CURVE SLOPE (PER RAD.) = 520624%000D OO 


MOMENTS OF TNERTIA(INZ Sees) 


it 


[T(X) Oe cOcoae(0 01 I(Y¥)) = 226436244) O2 


1(2) 26 VARs Leper ftXe) = Mieeeeassip Ol 


HOC mec. Geel Uiiese ore; ht 10> REFERENCE AXIS ORIGIN (IN) 


XBAR 4. (068655D-02 VORRe so aoe CO P635/0-05 


iH 


ZBAR 


Filet Lao 1D—O2 


Cn 
mh 





Anse 
‘ee ’ 





WING AND TAIL SURFACE PARAMETERS: 


~WING CHARACTERISTICS-—- 


ROOT CHORD (IN) = 621500000D 
TIP CHORD (CIN) = 1.3500000D 
WING SPAN (CIN) = 3.00000000 
WING AREA (IN2) = 1.2000000D 
LeE~ SWEEP (DEG) = 2-¢3500000D 
DIHEDRAL ANGLE (DcG) = 6.000000 0D 
LIFT CURVE SLOPE (PER RAD) = 425309000D 


“HORIZONTAL TAIL CHARACTERISTICS- 


ROOT CHORD (IN) = Se Df DOUU 0D 
TIP CHORD CIN) = 9. 3000000D- 
TAIL SPAN (CIN) = 1.3125000D 
TAIL AREA (IN2) = 3¢1100000D 
Reece vee? 1 UEG) = 22.59200000D 
DIHEDRAL ANGLE (DEG) = 0.0 

PhepecUmvE SEGRE (PER RAD) = 4202110000 
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00 
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00 
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~VERTICAL TAIL CHARACTERISTICS- 


ROOT CHORD (IN) = 6.5000000D 


TIP CHORD (IN) = 1.5000000D 


&6.375000CD 


TAIL HEIGHT (IN) 


TAIL AREA (IN2) = 2660000000 


BoEs  OWEEr BEG) 5 «00000 00D 


LIFT CURVE SLOPE (PER RAD) = 2-9891000D 
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LONGITUDINAL STABILITY DERIVATIVES: 


ge en aeae awe 2 ae ae ee OW ae ae ae SP ae ee 3S ee See ee ee eb ee ee ee ee 6 ee ee ee ae 


CM- = VARIATION OF PITCHING MOMENT COEFFICIENT 
Wot 
A = ANGLE OF ATTACK (CALFA) = PA Ch aRATE 


A(DOT) = VARIATION OF ALFA WITH RESPECT TO TIME 


CMQ = -1.49790/70D 01 CMA = -1.1440915D 00 


CMA(DOT) = -6,2308383D 00 CMA(TAIL) = -2.9448839D 00 


ELEMENTS CF STATE VARIABLE MATRIX A ARE: 


—-4.49702340-03 —2.4251559D-04 
1-Q000000D 00 0-0 


ELEMENTS OF THE B VECTOR ARE: 


-6.2438198D-04 
0.0 
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LATERAL STABILITY DERIVATIVES: 


_—_—=_e ae 8 ee ee ee eee ee eee eee ee ee ee ee 


CL- = VARIATION OF ROLLING MOMENT COEFFICIENT 
WITH - 
CN- = VARIATION OF YAWING MOMENT CCEFFICIENT 
WITH - 
= SIDESLIP ANGLE (BETA) 
= ROLL RATE R = YAW RATE 
DA = AILERON DEFLECTICN DE =" KUDDER “DEFEECTION 


lB == 1425372150-01 CNB = Te 5 1906501 0-02 
CLP = -4.3806752D-01 UNP = —2.5131716D-03 
CLK = 4.5663509D-02 CNR = ~-1,55841920-01 
CLDA = 6.9936720D-01 CNDA ="> 1 - Be / (OS2)-053 
CLDE =  5.5541463D-02 CNOF = -2.3746/359D-01 


ELEMENTS OF STATE VARIABLE MATRIX A ARE: 


-1.274/78840-01 1.1230146D-02 S52 ONG DSO 2 
-4.3405394D-03 -6..41 3000 60-03 5.05249 7207-03 C20 
0.0 =. VOCOOUGOD “O00 0.9 
1.0000000D 00 0.0 0.0 0.0 


ELEMENTS OF THE B MATRIX ARE: 


POS 4O( G99 7OL 1.38849610-02 
Sao oU4o020) O03 SLO et0s0 c 
0.0 0.0 
0.0 0.0 
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